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Abstract 

By applying algebraic techniques, we construct a two-parametric family of strictly 
isospectral Hydrogen-like potentials, as well as some of its one-parametric limits. An 
additional one-parametric almost isospectral family of Hydrogen-like potentials is 
also investigated. It is argued that the construction of a SUSY partner Hamiltonian 
using a factorization energy 5 less than the ground state energy of the departure 
Hamiltonian is unnecessarily restrictive. 
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There now exists a considerable amount of work relating to exactly solvable one- 
dimensional potentials in Quantum Mechanics (QM). The main interest has been to en- 
large the number of analitically solvable potentials using diverse techniques: the Darboux 
transformation |TJ, the Gelfand-Levitan formalism the standard and modified factor- 
izations [3-4], the supersymmetric Quantum Mechanics (SUSY QM) ||, etc. The un- 
derlying ideas of most of these procedures have been summarized in an algebraic scheme 
in which a first order differential operator intertwines two different Hamiltonians ||. In 
particular, the generation of one-parametric families of potentials isospectral to the tra- 
ditional ones [2-9] can be performed by means of this first order intertwinning technique. 
This method leads also in a natural way to the factorization of the implied Hamiltonians. 

Recently, a generalization of this technique, in which the intertwinning operator is of 
second order, has been used to derive a two-parametric family of potentials isospectral 
to the harmonic oscillator [10]. This second order intertwinning technique (SOIT) is a 
particular case with n = 2 of the n-order intertwinning technique which appears to be the 
best way to introduce the higher order SUSY Quantum Mechanics [11]. In principle, the 
latter makes possible to generate n-parametric families of Hamiltonians isospectral to a 
given Hamiltonian |12|, H| . 



In this paper we will generate a 2-parametric family of radial isospectral Hydrogen-like 
potentials by means of the SOIT. We shall show also that the one-parametric family of 
potentials derived by Fernandez j8| can be recovered from ours. As a final result, we will 
find a 1-parametric family of potentials having the same energy levels as the corresponding 
radial Hydrogen-like potentials except for the ground state energy level. 

The standard procedure to deal with Hydrogen-like potentials in QM reduces to solve 
the eigenproblem for a particle in a 1-dimensional effective potential Vi(r) = 1(1 + l)/r 2 — 
2/r, where I = 0, 1, 2, ... is the azimuthal quantum number and r is a dimensionless radial 
coordinate. By simplicity, instead of working with the standard radial wavefunctions R(r), 
we will work with the functions ip(r) = rR(r) with an inner product defined by (if), ip') = 
4-7T Jo +0 ° ip(r)xjj'(r)dr. As it is well known, the eigenvalues of the radial Hamiltonian Hi = 
—d 2 /dr 2 + Vi(r) (with fixed I) are given by: 

= £ lfc = --L-; k = 1,2,3,... (1) 

where I + k = n. 

We are looking for a Hamiltonian Hy = —d 2 /dr 2 + Vi>(r), such that the following 
intertwinning relationship is satisfied |J: 

He A = AHi, (2) 
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where A is a second order differential operator to be determined 



(3) 



Equations (2-3) and the explicit form of Hi and Hy lead to the following relations between 
Vi(r),V v (r),(3(r) and 7(7-): 



0' - y + (27(f) - /?' - fjj P 2 + 2c = 0, 



2 7 (r) = /? 2 - /?' - 2 



(4) 

(5) 
(6) 



where c and d are, in principle, arbitrary constants and the prime denotes derivative with 
respect to r. The key point becomes now to solve the non-linear second order differential 
equation (4) for f3(r). Let us first enforce that the operator A includes ai-\ai as a particular 
case, where a; = (—d/dr + l/r — l/l) is the standard factorization operator for the radial 
Hydrogen-like Hamiltonians 0. This condition permit us, without a loss of generality, to 
fix the constants c and d in (4-6) as 



(2/-1) 2 1 + (2Z-1) 2 

°~ 4/ 4 (/- l) 4 ' ~ 2Z 2 (/- l) 2 ' 

After including (7) in (4-6), it turns out that the general solution to (4) becomes: 



(7) 
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Z = 2,3,..., 



,r/I(J-l) 



1 - 



1(1-1) 



1 + 



(2/)! \l 
v 2 



x 2l e- 2x ' l dx 



21-1 



r x 2l e -2x/(l-l) 



dx 



(9) 
(10) 



(2Z — 1)! U-l/ Jo [Z(Z — 1) — ac] 
where v\ and 1/2 are integration constants. 

The asymptotic behaviour of (3(r) is given by (3(r) ~ (1 — 21) /l(l — 1), while it diverges 
as (21 — l)/r, when r — > + . This suggests to us to write the new potential Vi>(r) in (5) as 



V^ 2 (r) = V^ 2 (r) + 2a'(r), 1 = 2,3,..., 



where V = I — 2 and a(r) = /3(r) + (1 — 21) /r, is an appropriate function that makes 
evident the limit Vi^ 2 (r) — > VJ„2( r )> when r — ► +00 and r — > 0. The parameter domain 
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for which a'(r) is free of singularities is given by vi, v 2 G (—00, 1); inside this parameter 
region, the new two-parametric family of potentials VJ_ 2 (r) has the same singularities as 
Vz_ 2 (r). Furthermore, for v\ = v 2 = we have a' — and Vi- 2 (r) = VJ_ 2 (r). This means 
that Vi^ 2 (r) not only governs the asymptotic behaviour of V^_ 2 (t) but it is a member of 
the family (11). 

Now, from equation (2) it becomes apparent that the operator A acting on the eigen- 
functions {ip n i( r )} of Hi provides eigenfunctions ■0 r j ) ;_ 2 (r) oc Atp n i(r), I = 2,3, ... of H_ 2 
with eigenvalues E ik = E [k , i.e.: 

fe- 2 (r) = ^-1K A ^ (r)- (12) 

/ ( n 2_/2)( n 2_ /2 + 2 /_ 1) 



However, the set {ip n ,i-2ij')} is n °t yet complete in the Hilbert space of square- 
integrable functions 7i. This is clear if we try to answer the following questions: are 
there functions ^_ 2 , e (r) orthogonal to all the -^ ni ;_ 2 (r)? If so, are they eigenfunctions of 
Hi_ 2 l In order to answer these questions, let us assume that the set {^_ 2 , e } exists, then 

(^_ 2 , e (r),Vv_ 2 (r)) « $i-2,e( r l A M r )) = (A^i- 2 , t {r),^ nl {r)) = 0. 

Due to the fact that {Vvu( r )} is a complete set in H, the kernel of the second order 
differential operator A' is a 2-dimensional subspace 7i t orthogonal to all the ■0 n ^_ 2 (r), 
/ = 2,3, ... Let us write them as ^_ 2je (r) = c exp[J where c is a constant and 

/(x) is to be determined. The equation A^_ 2je (r) = can be rewritten as: 

f (r) - P(r)f(r) + f{r) - /3'(r) + 7 (r) = 0. (13) 
This Riccati type differential equation has a general solution given by 

/(r) = y - l - + /3(r) + ^ ln[c l9l (r) + c 2 <? 2 (r)], (14) 

where <?i(r) and <? 2 (r) are given in (9-10), and ci and c 2 are constants. The generic kernel 
element is given by ^- 2 , e (r) = Co^_ 2;0 (r) + C_i^_ 2 _i(r), where 



and 

^-2,-i(r) 
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are both eigenfunctions of H1-2 with eigenvalues £j-2,o = — 1/Z 2 , ^-2,-1 = _ ~~ 1) 2 > 
/=2,3,... respectively. In (15-16), W{g\,g<i) = g[(r)g2(r) — g\{r)g' 2 {r) represents the 
Wronskian of g\ and g^- 

Notice that {$1-2 -1, $1-2,0, $>n,i-2, 1 — 2,3, ...} is now a complete set in TL, and their 
elements are eigenfunctions of Hi_ 2 , I = 2, 3, with eigenvalues: 

111 

El - 2 > k = ~p' ~WW ] fc = 1 ' 2 ' 3 '- (17) 

Comparing with (1), one can see that this spectrum is identical to that of H1-2, 
i.e., Ei_2,k = Ei-2,k- Hence, the Hamiltonian if;_2 is strictly isospectral to -£^_2, and 
because VJ_2(r) depends of two free parameters v\ and V2, a new two-parametric family of 
isospectral Hydrogen-like potentials has been generated. Some particular cases are worth 
discussing in more detail. 

Firstly, for v\ — v 2 = 0, the well known Hydrogen-like potential Vi-2^) is always 
recovered because it is a member of the family V^_2( r )- 

Let us now to take z/i = iv2 = 0), in this case the 2-parametric family VJ_2(r) is a new 
1-parametric family which becomes to a Hydrogen-like potential when z/2 — > iv\ — > 0). 
In particular, let ui — 0, and 



^2 



'Z — 1 \ ^ 1 (2/ 



v 2 / 7i-i 
In this case, the function a(r) becomes 

1-2/ r 2/-2 e -2r/q-l) 

a(r) ~ Z(Z-l) + 7z _ 1 -/ r x 2 '- 2 e- 2 ^- 1 )rfx' / ~ 2 ' 3 '- 

If we change in (11) I by Z + 1 with a(r) as it is defined above, we get the following 
one-parametric family of potentials: 

^ lW ^_ l(r ) + 2 _|___j, i = 1 , 2 ,... (18 ) 

which was generated by Fernandez in 1984 j8|. Notice that, in the case when 7/ — > 1/4, 
the family (18) gives the particular case derived by Abraham and Moses |7]]. Moreover, 
when 7; — > 00 (z/ 2 — * 0), we have V^_i(r) — > VJ_i(r), just as we have proposed. 

The previous potentials can be seen as deformations of V/_i(r) induced by the second 
term in (18), which does not change the behaviour of V5_i(r) at the ends of the interval 
[0, 00), but can produce important modifications inside. In particular, there is the possi- 
bility of creating one additional well in VJ_i(r) whose depth and position can be changed 
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by varying 7/. On the other hand, the two-parametric family of potentials (11) admits 
also the previous interpretation, but more freedom is given by the deforming term (it 
depends on two parameters instead of just one as in the previous case). Hence, we get 
the possibility of introducing now two wells, one of them with its minimum placed around 
the global minimum of Vz_ 2 (r) and the other one pushed further out. As in the previous 
one-parametric case, the depths and positions of the two wells can be modified by chang- 
ing the two parameters v 1 and v 2 . This is illustrated in Figure 1, where two members 
of the family Vi(r) of (11) are depicted, together with the undeformed potential Vi(r) 
(dashed line). We have plotted also in Figure 2 the corresponding probability densities 
for the two energy levels E1-1 = —1/4, E lfi = —1/9 with v\ — v 2 — —10. For the lowest 
level the probability has a maximum around the left well. On the other hand, the first 
excited state has two maxima, the highest one centered around the right well, while the 
lowest one is situated around the left well. When we go over the higher excited states, 
the probabilities resemble more and more the corresponding Hydrogen-like densities. 

Up to now, departing from the Hydrogen-like potential VJ(r) we have generated a two- 
parametric family of solvable potentials Vi- 2 (r), with the same spectrum and singularities 
as Vi-iif). Now, as the intertwinning operator A is of second order, it is interesting to 
look for its possible factorization in terms of two first order differential operators b±, and 
b 2 , i.e.: 

A = b 2 h, b j = -^- + w j (r), j = 1,2. (19) 

This leads to 7(7") = w' 2 + w\w 2 , and wi(r) = (5{r) — w 2 (r), where (5{r) is given by (8), 
and w 2 has the form given in (14). Hence, W\(r) takes the form: 

w±(r) = - - y - ^- ln[ci0i(r) + c 2 g 2 (r)\. (20) 
r I ar 

It is now clear that when solving (13), we simultaneously have gotten the solutions to 
the equation A^i_ 2i6 = 0, as well as the factorizations of the operator A. There is a 
continuous family of factorizations because when we change the values of c\ and c 2 in 
(20), we are simultaneously changing the operators b\ and b 2 , but maintaining fixed their 
product A = b 2 b\. 

A first consequence of this factorization arises after rewritting equation (2) as H\^ 2 b 2 b\ = 
b 2 b\Hi. Suppose now that there is a Hamiltonian H*_ x = —d 2 /dr 2 + VJ* x such that 
#,*_!&! = hHi. Thus & 2 #f_i&i = hhHi, and we get H X - 2 b 2 = b 2 HJ f _ 1 . Therefore, Hf_ Y 
could be considered as an intermediate Hamiltonian between Hi, and iJz_ 2 . Hence, the 
SOIT can be seen as the iteration of two first order intertwinning transformations. 
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We notice here that from the very beginning we are labeling with the subindex / — 1 the 
intermediate Hamiltonian H*. This is a consequence of the further calculations leading 
to a centrifugal term for V* with exactly that index. 

In order to ensure that the first order intertwinning relationship H^bi = b x Hi would 
be satisfied, the functions W\ and V*_\ must satisfy some restrictions. The key one becomes 
the following Riccati equation: 

-w' 1 + (w 1 ) 2 -V l + 5 1 = 0, (21) 

complemented with the typical SUSY relationship V l *_ l = V\ + 2w[, where Si is a constant 
to be determined. We notice also that a first order intertwinning relationship of the kind 
H*_ibi = biHi leads in a natural way to the factorization of the Hamiltonians Hi and 
Hl_ L : Hi = b\bi + S u and H*_ 4 = bib[ +Si. By a similar argument, the Hamiltonians H*^ 
and Hi_2 become factorized in terms of b 2 and b\: H* = b 2 b 2 + 82, and H = b 2 b 2 + S 2 . 

In order to determine the intermediate Hamiltonian Hf_ x , we must find inside the 
family (20) a member obeying equation (21). In fact, for ci = and c 2 = 1, we have 
Si = — l) 2 , / = 2,3, .... Thus, the potential V i I 1 (r) can be written as: 



W2Y ~ 9292 



(22) 



{92? 

The parameter domain for which the family (22) has the same singularity as VJ_i(r) is 
given by v 2 G (l,oo). The eigenf unctions of H*_±, I = 2,3,... are given by -0*— 1 — 1 cx 
r l e~ r / l /g 2 , and ip*^ oc biip nh with eigenvalues E?^ ^ = -1/(1 - l) 2 , and Ef_ 1>k = E lk , 
k = 1,2,..., respectively. Notice the unusual absence of the state corresponding to E^ = 
-l// 2 . A direct comparison of the spectra shows that Vj*_i{t) is almost isospectral to V1—1, 
the difference resting on the ground state level position. 

The next first order intertwinning transformation gives a different factorization of H^_ l 
and some interesting new results. The absent energy level ^-2,0 — is added now 
to the spectrum of H*_ x in order to generate Hi_ 2 . But this means that the factorization 
energy in this second step is greater than the ground state energy level of H^_ 1 , and this 
naturally fills the hole generated by the first factorization. The eigenfunctions of Hi_ 2 , 
in terms of those of H*_ u I = 2,3, are given by {^-2,-1 oc b 2 ip*_i_i : ^-2,0, VW-2 oc 
b2^n 1-1}, where ^-2,0 oc exp(— / w 2 (x)dx) is the eigenfunction associated to the 'missing' 
energy level E?_ 1 = Ei_ 2 . A direct calculation shows that this set of eigenfunctions is 
the same as that derived by means of the SOIT. 

Concluding Remarks. In this paper we have shown that the Second Order Inter- 
twinning Technique allows one to derive a two-parametric family of isospectral Hydrogen- 
like potentials. The iteration of two first order intertwinning transformations leads to the 
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same results but gives additional information. Thus, against the standard statement of 
SUSY QM, there are cases where a factorization energy greater than the ground state 
energy of the departure Hamiltonian leads to a physically acceptable SUSY partner. A 
deeper discussion of these first order intertwinning cases will be given elsewhere [ I3| . Fi- 
nally, all the potentials presented here have the same kind of singularity at a fixed point 
(r = 0) as the initial Hydrogen-like Hamiltonian. The case when the intertwined poten- 
tials are non-singular has been illustrated by Fernandez et. al. for the harmonic oscillator 



10, [12| , while the case when the SUSY partner of the oscillator potential has a mov- 



able singularity has been successfully interpreted [ 14[ . The corresponding problem for a 
family of isospectral Hydrogen-like potentials with a different singularity as those of V/(r) 
in r = is open. 
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Figure 1: The Hydrogen- like potential VJ_2(V) (dashed line) and two members of the 
family Vi^ij') (solid line) with 1 = 3. Here we have chosen equal values for v\ and z/ 2 , but 
it is not a restriction (see eq 11 and bellow). Notice that in the limit case when z/j — > 0, 
i = 1,2, the solid lines resemble more and more the dashed one. 
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Figure 2: Behaviour of the probability densities of the two first energy levels -1/4 and 
-1/9 of H 1: with vi = v 2 = —10. The levels are indicated by the dashed lines 
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